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We use gauge/string duality to investigate the free energy of two static color sources (a heavy 
quark-antiquark pair) in a Yang-Mills theory in strongly interacting matter, varying temperature 
and chemical potential. The dual space geometry is Anti-de Sitter with a charged black-hole to 
describe finite temperature and density in the boundary theory, and we also include a background 
warp factor to generate confinement. The resulting deconfinement line in the \i — T plane is similar 
to the one obtained by lattice and effective models of QCD. 
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Consider QCD in a four dimensional Euclidean space- 
time and in nuclear matter, and two static color sources, 
an infinitely heavy quark and an infinitely heavy anti- 
quark, at distance r from each other. It is interesting to 
investigate how the free energy of such a system behaves 
against variations of temperature and chemical potential. 

The study of a strongly coupled Yang-Mills theory, 
such as QCD, is a challenge in spite of the methods de- 
veloped so far to deal with it (lattice simulations, mod- 
els and effective field theories). The formulation of the 
gauge/gravity (or Anti-de Sitter/conformal field theory) 
correspondence [IH1] has suggested to face this problem 
through the identification of a suitable higher dimen- 
sional gravity dual. Due to the strong/weak nature of 
the duality, the gravity dual is a weakly coupled theory 
defined in a higher dimensional curved spacetime; since 
QCD is not conformal, a mechanism for breaking confor- 
mal invariance must be included in the dual model. 

This holographic framework can be adopted to ana- 
lyze finite temperature and density effects. In a bottom- 
up approach, we use the soft wall AdS /QCD model, a 
five dimensional model formulated on AdS§ spacetime, in 
which linear confinement at zero temperature and den- 
sity is obtained by inserting a background warp factor 
[5HZ], bringing a mass scale related to A QCD . To describe 
the boundary theory at finite temperature, a black-hole 
is included in the five dimensional space, whose horizon 
position represents the (inverse) temperature 0]. 

On the other hand, in QCD the effect of finite 
quark density is introduced by adding the term Jo — 
(x)ip(x) to the Lagrangian in the generating func- 
tional, so that the chemical potential fi appears as the 
source of the quark density operator. According to the 
AdS /CFT correspondence, the source of a QCD opera- 
tor in the generating functional is the boundary value of 
a dual field in the bulk; therefore, the chemical potential 
can be considered as the boundary value of the time com- 



ponent of a (7(1) gauge field Am dual to the vector quark 
current. Under the ansatz A = A a (z) (z is the fifth 
holographic coordinate) and Ai = A z = (i = 1,2,3), 
one can find a solution of the equations of motion of a 
5D gravity action with negative cosmological constant 
interacting with an electromagnetic field: the solution is 
known as the AdSVReissner-Nordstrorn black- hole met- 
ric, and describes a charged black hole interacting with 
the electromagnetic field [8j [9] . We use this (Euclidean) 
metric in our model: 

ds 2 = - 2 f f(z)dt 2 + dx 2 + jj-r\ , 

M = l-^+ q 2 z^z' + q 2 z\ (1) 

where R is the AdS radius, q the black-hole charge and 
Zh the position of the horizon, defined by the condition 
f( z h) — 0. The e c z term, characterizing the soft wall 
model, distorts the metric and brings the mass scale c 
HO] . The positive sign in the exponent is chosen, following 
ref.[B], to obtain a confining QQ potential at T = 0. 
The profile introduced in [7] , with the negative sign, also 
produces linear Regge trajectories and bound states for 
light hadrons, but does not provides an area law for the 
Wilson loop, as discussed, in the framework of light front 
holography, in [11 . The boundary z = represents 
the four dimensional spacetime on which the Yang-Mills 
theory is defined. 

The Hawking temperature of the black-hole is 

47r dz z=z h **h V 2 J 



Issues related to the choice of the dilaton profile are discussed in 



2 



where Q = qz\ and < Q < moreover, the Eu- 
clidean time is a periodic variable, with period f3 = 1/T 
[H H3] . Studies of this kind of models for q = can be 
found in [13] . 

The relation between the parameters in the metric ([I]) 
and the chemical potential can be obtained observing 
that, on dimensional grounds, the low z behavior of the 
bulk gauge field A (z) is 



A Q (z) =i-i- 1]Z 2 



(3) 



with rj — nq and n a dimensionless parameter in our 
model. Together with the condition that Aq vanishes at 
the horizon, A (zi l ) = 0, Eq.(|3]) determines a relation 
between \x and the black-hole charge: 



[i = Kqz h 



(4) 



The parameter c in the warp factor in sets the scale: 
in the following, we analyze quantities expressed in units 
of c, and comment on the numerical results at the end. 

Let us now turn to the problem formulated at the be- 
ginning of this study. At finite temperature, the free en- 
ergy F(r, T) of an infinitely heavy quark-antiquark pair 
at distance r can be obtained in QCD from the correla- 
tion function of two Polyakov loops: 



(P(fi)P t (f 2 )) = e"* F(r ' T)+7(T) 



(5) 



with r — \xi — x 2 | and j(T) a normalization constant. 
Moreover, the expectation value of a single Polyakov loop 



(T) 



(6) 



(F°°(T) = F(r = 00, T) and neglecting the normal- 
ization) is the order parameter for the deconfinement 
transition in a pure SU(N) theory [15] . Within the 
gauge/string duality approach, we can attempt a calcu- 
lation of the expectation values in ^ and ^ consider- 
ing string configurations in the 5D manifold having the 
Polyakov loops as boundary, looking at the configurations 
of minimal surfaces and computing their worldsheet ac- 
tion. The worldsheet action is the Nambu-Goto one 

d 2 ^det [g MN (d a X M )(d b X N )} (7) 



2-na' 



for a string with endpoints attached to the positions of 
the heavy quark and antiquark, x = ±r/2 at z — 0, so 
that QB H2 



F(r,T) =TS KG . 



(8) 



X M, s are the coordinates of the five dimensional space- 
time, gMN the metric tensor associated to the line el- 
ement 0, and the determinant is taken over the in- 
dices a, b = 0,1 labeling the worldsheet coordinates £°. 
We parametrize them with £° = t and = x, and 
look for a static solution z(x) satisfying the conditions 
z(x = 0) = z , z'(x = 0) = and z(x = ±r/2) = 0, the 
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(a) confined 



(b) deconfincd 





(c) confined 



(d) deconfincd 



Figure 1: String configurations corresponding to the confined 
(a,c) and deconfined phase (b,d). In (a,c) there is one string 
with endpoints at the positions of the static color sources, 
x = ±r/2 at z — 0, which does not intersect the horizon 
(dark cylinder in (c)); in (b,d) there are two strings stretched 
between z = and the horizon z — Zh, with x(z) = ±r/2. 



prime denoting the derivative with respect to x. Exam- 
ples of such configurations are shown in Fig. l(a)| 1 (d) 



From ([7]) and ^ we obtain an expression for the free 
energy: 



F(r,T) 



dx- 



r/2 



\fm + {z'f 



(9) 



with g = -^r . Moreover, an equation of motion follows 
from ([9]), with the first integral 



H = 



m 
\fm + {z'f 



(10) 



which allows us to express F(r, T) in terms of zq and 
/o = f(zo). Defining v = z/zq, after having subtracted 
the UV (v — > 0) divergence (corresponding to subtract 
the infinite quark and antiquark mass in four dimensional 
QCD [IB]), we obtain 



F(X) 




(11) 



where F = F/c, A = cz and 
r(v) = aI 



fo 



v 4 e 2\ 2 (l-v 2 ) 



(12) 



The distance f — or can also be expressed parametrically, 



f(A) = 2\Jh 



dv 



v 2 e A ) 



r(v)f(z v) 



(13) 
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Figure 2: Interquark distance f(A) in the confined (a) and 
deconfined phase (b). In (a) (here for p = 0.5 and T = 
0.82T*), for ^ A ^ Ai, f(A) spans the whole real axis (red 
solid line); as A2 ^ A ^ zh, r(X) monotonically decreases from 
a finite value r(\2) to r(zu) = (blue dashed line). The free 
energies corresponding to the two branches are plotted in (c). 
In (b) (here for p — 0.5 and T = 1.65T*) f (A) has a maximum 
for Ai = A2 = A*; the free energy is plotted in (d). 



so that, from (111 and (13) we can write the free energy 
as a function of the temperature T = T/c, for a range 
of values of the chemical potential p — fi/c. We fix the 
parameters k and g to one, postponing to the end the 
discussion on the numerical results. 

As for p — [17], also for finite chemical potential 
the distance r(A), written parametrically in ( 13 ), has two 



branches, shown in Fig. 2(a) For ^ A ^ Ai, f increases 
monotonically from f (0) = 0, spans the whole real axis 
and diverges when A — > Ai. For A 2 ^ A ^ Zh = cz^, f 
monotonically decreases from a finite value r(Aa) down 
to f(zh) — 0; for Ai < A < A2, f has an imaginary part. 
Selecting the branch ^ A ^ Ai, we can comp ute F for 
all distances, with the result depicted in Fig. 2(c) As 



shown in the same figure, the values of the free energy 
corresponding to the branch A2 ^5 A Zh are larger than 
those for Q ^ A ^ Ai. The resulting free energy can be 
parameterized similarly to the Cornell expression for the 
static quark-antiquark potential: 



F = -- +b + a 2 f 
r 



(14) 



The coefficient a turns out to be essentially independent 
of T and fi; on the other hand, at fixed T = 0.05, the 
parameter b increases, from b = —0.073 at p = to 
b = —0.051 at pL = 0.6, while a decreases from a = 0.64 
to a — 0.62 when p is changed from p = to p = 0.6. 

For fixed values of the chemical potential p = p* , in- 
creasing T, the two points Ai and A2 along the A axis 
get closer to each other, and at a certain value T* they 
coincide: Ai = A 2 = A*. For f > f*, f(A) is real for 




Figure 3: Free energy in the confined phase, p = 0.5, f = 
0.82T* (blue solid line), and in the deconfined one, for p — 0.5, 
f = 1.23T* (green dot-dashed) and p = 0.5, f = 1.65T* (red 
dashed). For p = 0.5 the critical temperature is T* = 0.122. 



^ A ^ Zh, a nd never exceeds the value f (A*), as shown 
in Fig. 2(b) at the point (p*,T*) in the p — f plane 
f(A) is always bounded. Therefore, taking p — p * anc ' 
T ^ T* , F can be evaluated only for f ^ r(A*), Fig. 2(d) 



For r > f(A*) there is another string configuration, i.e. 
the one with two strings stretched between the bound- 
ary z = and the black- hole horizon z = Zh, plotted in 
Figs. 1(b) and l(d)| 18J. This configuration holographi- 
cally represents two deconfined quarks. To evaluate the 
corresponding action we choose a different parametriza- 
tion of the worldsheet: £° = t and £ x = z, obtaining the 
regularized free energy 



1 

Zh 



Zh dz 



(e--!) 



+ C(P,f) (15) 



but 



with z — cz; £(p,T) is a constant related to the regular- 
ization procedure, and fixed by matching F(f*) = F° 
In this case the free energy does not depend on f, 
only on the temperature and chemical potential. 

As a result, the free energy is depicted in Fig. [3| At 
large values of f, for a given value of the chemical poten- 
tial p there is a temperature T*(p) such that for T less 
than T* the free energy linearly grows with f, while for 
T ^ T* the solution coming from the connected string 
configuration (Fig. 1(a)) is only valid up to f = r(A*); 



from that point on, F is constant (15). This behavior of 



the free energy with distance and temperature has been 
observed in quenched lattice simulations at p — |19j . 

In the p — T plane, the curve defined by the points 
(ft, f* (p)) is plotted in Fig.[i] This can be considered as a 
picture of the deconfinement transition in the QCD phase 
diagram. For low values of T and \i there is a confined 
phase characterized by a linearly increasing interquark 
potential. Outside this region, a finite amount of energy 
is sufficient to separate the quark-antiquark pair, as ex- 
pected in a deconfined phase. The picture agrees with 
the diagram obtained by, e.g., Nambu-Jona-Lasinio |20j 
and other effective models [2Tj . 

In the same holographic framework, the Polyakov loop 
(V) can be computed by Eq.([6]). In the confined phase, 
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Figure 4: Deconfinement transition line in the pb — T plane. 
The line divides the plane in two regions, a hadron phase near 
the origin, and a deconfined phase beyond the curve. 



F(f) diverges as f — > oo, so (P) vanishes, while in the 
deconfined one it is determined by F°° from the config- 

The result is depicted in 



uration in Figs. 1(b) and 1(d) 
Fig. [5| for each value of fi, the Polyakov loop vanishes as 
T < T*; it starts growing at f = f* [22]. The smooth 
increase can be interpreted as a continuous transition be- 
tween the hadron and the deconfined phase. 

As pointed out in [9], the parameter k scales as k ~ 
y/ N c and its numerical value is model dependent. If we 
take the limit N c — > oo, the deconfinement line in Fig. [4] 
becomes flat, with T* not depending on /I, in agreement 
with the result found in |23j . Considering a finite k in 
Q, the deconfinement line has the shape shown in Fig. [4] 



In order to discuss numerical results, we need to choose 
the values of the parameter k in Q and of the mass scale 
c appearing in the warp factor. We fix k « 1/2 from the 
ratio f c /p, c « 0.5 found in 20J, with f c = f* at £ = 
and (l c = fi* at T = (k could also be determined in this 
soft wall model studying different observables, as done in 
[9] for the hard wall model). The mass scale c induced 
by the warp factor can be fixed from the result T c ss 170 
MeV gl] : c « 850 MeV. On the other hand, the value 
of this scale from the p meson mass is c « 670 MeV [5], 
which gives T c w 134 MeV and (j, c w 248 MeV. 2 

The numerical values of the critical chemical potential 
are a remarkable result of our study. The other result 
is that the gauge/string duality approach can be used 
to obtain information on the deconfinement transition in 
QCD, finding an agreement with the outcome of effec- 
tive theories (even if the model is not sensitive to the 
structure of the vacuum) . The success of the holographic 
framework is noticeable, considering the difficulty of lat- 
tice QCD to explore the phase diagram along the axis 
of the chemical potential, and opens interesting perspec- 
tives for the analysis, for example, of the superconducting 
phases of QCD. 



CP) 



0.2 



Figure 5: Polyakov loop (V) versus t (normalized at its value 
for T — > oo) for p, — 0.2 (blue, dashed, right curve) and 
ft — 0.7 (red, dotted, left curve). 
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